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We study the electronic spectral properties in two examples of strongly interacting systems: a 
Mott-Hubbard insulator with additional electron-boson interactions, and a polaronic semiconductor. 
An approximate unified framework is developed for the high energy part of the spectrum, in which 
the electrons move in a random field determined by the interplay between magnetic and bosonic 
fluctuations. When the boson under consideration is a lattice vibration, the resulting isotope effect 
on the spectral properties is similar in both cases, being strongly temperature and energy dependent, 
in qualitative agreement with recent photoemission experiments in the cuprates. 

PACS numbers: 71.38.-k, 79.60.-i, 74.72.-h 



I. INTRODUCTION 

A general feature of strongly interacting electron sys- 
tems is a sizeable suppression of the metallic character, 
accompanied by a redistribution of the spectral weight 
from the region close to the Fermi level towards higher 
energies. In the presence of a strong electron-electron 
repulsion, for example, incoherent excitations arise far 
from the Fermi energy in the so-called upper and lower 
Hubbard bands, at an energy scale which is ruled by the 
strength U of the interaction. A similar behavior is also 
found in systems with strong electron-lattice coupling, 
where the spectral weight is transferred to a broad peak 
located around the polaron binding energy Ep. 

Recent photoemission experiments in the high-Tc 
cuprates 0, 0, 0, H, which are generally described in 
terms of purely electronic models due to the proxim- 
ity to a Mott insulating phase, have revealed the ex- 
istence of an important electron-lattice coupling. This 
has given rise to intense theoretical work focusing on the 
excitation spectra in the presence of electron-boson in- 
teractions, both in models with 0, 0, 0] and without 
pi. IsL fToL [ill IT2I IT^ I electronic correlations. Of particular 
interest are those works which focus on isotope effects 
(IE), since they can disentangle the properties which are 
directly related to the coupling to the lattice degrees of 
freedom. For example, it has been found that an anoma- 
lous isotope effect on the effective mass arises at the po- 
laron crossover, signaling the breakdown of the Migdal- 
Eliashberg adiabatic approximation ^4|. An enhance- 
ment of the IE is also expected in the proximity of a 
Mott metal-insulator transition, i.e. at intermediate val- 
ues of the electron-electron repulsion 1.5$ . On the other 
hand, a strong electronic repulsion suppresses the effects 
of the lattice dynamics on the low energy excitations 
near the Fermi level. These persist at high energy, where 
they amount essentially to a broadening of the electronic 



spectra Q- 

In this work, we provide an approximate theoretical 
framework based on the Coherent Potential Approxima- 
tion (CPA), which qualitatively describes the high-energy 
excitation spectra resulting from both electron-electron 
and electron-boson interactions, treating the fluctuations 
of magnetic and bosonic origin as a local static disorder. 
Although simplified, the present analytical treatment is 
able to account for the existence of broad incoherent 
peaks at high energy, which are a common characteristic 
of strongly interacting electronic systems. 

We analyze two extreme cases which can be impor- 
tant for our general understanding of the problem: (i) 
a Mott-Hubbard insulator at half-filling, in the presence 
of an additional local interaction of the electrons with 
dispersionless bosons and (ii) a polaronic semiconductor, 
where the physics is solely determined by the electron- 
boson coupling. In the former case, the validity of our ap- 
proach relies on the separation of electronic and bosonic 
energy scales, which is achieved due to a strong electron- 
electron repulsion. In the latter case, the present CPA 
results are controlled by direct comparison with an ex- 
act solution obtained by Dynamical Mean Field Theory 
(DMFT) 0. 

The consequences of the electron-boson coupling on 
the dispersion and width of the high energy features in 
the spectral function are calculated in both situations, 
in section II and III respectively. Particular attention is 
devoted to the effect of a shift of the boson frequency, 
as can be achieved through an isotopic substitution if 
the bosonic mode that couples to the electrons is a lat- 
tice vibration, or more generally if it is a collective mode 
with a sizeable lattice component (in which case a non- 
trivial IE can still arise, provided that the frequency of 
the boson is modified by the isotopic substitution) . In the 
present approximation, the IE on the high energy part 
of the spectrum turns out to be qualitatively similar in 
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the Mott-Hubbard insulator and in the polaronic semi- 
conductor, being mainly determined by the strength of 
the electron-boson coupling (albeit slightly enhanced by 
the presence of electronic correlations): in both cases it 
is strongly temperature and energy dependent, in agree- 
ment with the recently measured IE in the hi gh t em- 
perature superconductor Bi 2 Sr2CaCu208+,5 [TH Il8|. A 
tentative analysis of the experimental results, performed 
in section IV, is compatible with the existence of a mod- 
erate electron-phonon coupling in the cuprates. 



II. ELECTRON-BOSON COUPLING IN A 
MOTT-HUBBARD INSULATOR 

A. Coherent Potential Approximation 

We study the following Holstein-Hubbard Hamilto- 
nian: 
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where electrons in a band with dispersion e k mutually 
interact through an on-site repulsion U, and are coupled 
locally to a dispersionless bosonic mode of frequency too, 
with a strength g. We shall set the energy units such 
that h = ks = 1- We calculate the spectral properties 
of the above model at half filling, in the framework of 
the Coherent Potential Approximation (CPA) [l9l |20| . 
This is suitable for the Mott insulating phase at large 
U, but also gives a fair description of the high energy 
incoherent excitations in the correlated metal, provided 
that the upper and the lower Hubbard bands are well 
separated from the low energy quasi-particle peak |2l| . 
The latter, however, is not accessible within this theory, 
and the position of the chemical potential remains unde- 
termined except at half-filling. Additional features that 
are not included in the present description are the coher- 
ent excited states that arise at low temperature, (even 
in the insulating phase), due to the quantum nature of 
the magnetic excitations, in an ener gy rang e J ~ t 2 /U 
around the Hubbard band edges [221 l23t |24| . and whose 
dispersion is srongly reminiscent of spin density waves. 

Bearing these limitations in mind, the success of the 
CPA is that, despite its formal simplicity, it correctly 
accounts for the high-energy scattering by the ran- 
domly distributed magnetic moments. The momentum- 
integrated Green's function has the form 20]: 
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(1) 



where Go is an effective propagator which takes hopping 
processes into account. It can be eliminated by introduc- 



ing a local self-energy £(u;) through the following self- 
consistency condition: 



G( W ) = V 
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(2) 



(in the "atomic" limit, Go = 1/w and the usual single 
site propagator is recovered 20]). For each frequency uj, 
we are left with a system of two equations for the two 
complex unknowns G(ui) and E(w). 

In the presence of a local electron-boson interaction, 
equation Q can be generalized by introducing an addi- 
tional field y, which accounts for the random distribution 
of electronic energies due to the fluctuations of the bosons 
at different sites 



GM = J ^P{y) 
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In a system where the on-site electron-electron repul- 
sion directly competes with the attraction induced by 
the bosons, the formation of bipolarons (and the result- 
ing strong anharmonicities in the boson field) [2^, |2^| is 
prevented provided that U is much larger than the po- 
laron binding energy Ep = g 2 /coo ^Jllj}- In this case, 
we can assume that the boson field obeys the following 
gaussian distribution, 
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The variance 



u = EpujQ coth 
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(4) 



(5) 



is determined either by the thermal fluctuations (at high 
temperature, T 3> wq) or by the quantum fluctuations of 
the bosons (at low temperature, T <§; u>q). 

Let us emphasize that we have neglected boson renor- 
malization effects which, even in the presence of a strong 
electron-boson interaction, are suppressed by a suffi- 
ciently strong repulsion U 0, |2^] . Such effects can in 
principle be included by taking ujq in eq. (J3J) equal to the 
measured boson frequency rather than the bare parame- 
ter of eq. Q). However, the present scheme is expected 
to break down in situations where the bosonic and elec- 
tronic energy scales become comparable (i.e. Ep ~ U), 
in which case the properties of the bosonic field depend 
crucially on the electron correlations. 

Equations © and © form the backbone of our theory. 
The latter can be further explicited by evaluating the 
Hilbcrt-transform of the gaussian 



dyP(y)- 
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where W is the complex error function [28|]. Note also 
that, since the self-energy E does not depend on mo- 
mentum, the details of the band dispersion efc enter in 
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equation J5J only through the corresponding density of 
states. For the sake of simplicity, we shall consider a 
band dispersion with a semicircular density of states 
of half-width D, for which the self-consistency equation 
reduces to Gq 1 = lu — D 2 G/4. However, since they 
rely on a momentum independent quantity (the local self- 
energy), the results are quite independent on the choice 
of the band dispersion. 

The spectral function, which is the quantity of interest 
in the present work, is defined as 
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Its momentum integral 
N*(u) = -iJmG(w). 
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B. Mott-Hubbard insulator 
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FIG. 1: (color online) A sketch of the spectral density (the 
momentum-integrated spectral function). Top: in a Mott- 
Hubbard insulator, for U/D = 3 without (dashed green) 
and with (dotted blue) electron-boson interaction (Ep/D = 
0.9, u>q/D = 0.1). Bottom: in a polaronic system (U = 0) 
with the same electron-boson parameters. The green dashed 
curve here is the noninteracting band, the black solid line 
is the exact DMFT result [T^l . and the blue dotted curve 
is the approximate result based on equation 11311 . Vertical 
lines mark the reference energy of incoming particles. In both 
panels the shaded area represents the low energy part of the 
spectrum. 

For large U, and in the absence of electron-boson inter- 
actions, the set of equations and J2Jl can be rewritten 
as 
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The spectral density N*(u>) ~ D 2 /2 - {uj ± U/2) 2 

is illustrated in figure l.a, and consists of two bands 



of reduced width ~ 2D/^/2 separated by a gap A ~ 
U/2 — D/^/2 (the + and — sign are for the lower and up- 
per Hubbard band respectively). Such bands represent 
incoherent states that are strongly scattered by the dis- 
ordered magnetic moments. The corresponding spectral 
function A(ek,oj) exhibits broad peaks, whose width is 
comparable with the bandwidth itself: it is determined 
by the scattering rate T(u) = —ImY>(uj) ~ D 2 N*(uj) 1 
which is proportional to the spectral density, and is there- 
fore strongly energy dependent. As a consequence, the 
peaks sharpen (and become asymmetric) when approach- 
ing the band edges, as can be seen in the energy scans at 
constant momentum, in figure 2. a (dashed green curve). 
Their dispersion, defined as the locus of the maxima of 
A(ek,u)) at constant uj (the so-called momentum distri- 
bution curves, MDC) can be obtained by solving the 
equation Ek — ReT,(Ek) — and is roughly given by 
E k = -U/2 + e fc /V2 (see fig. 2.b). 

It should be noted that the slope of the dispersion of 
such high energy features deviates significantly from the 
unrenormalized Fermi velocity vf- For this reason, in 
strongly correlated electron systems, some care should be 
taken when extracting the function ReT,(uj) from experi- 
mental data under the assumption that vf tends asymp- 
totically to its unrenormalized value at sufficiently hig h 
energy, as is customary in the field (see also ref. |30|h 



C. Electron-boson coupling 

The coupling to a bosonic mode leads to several mod- 
ifications of the picture described above. First of all, the 
fluctuations of the site energies lead, through the relation 
Q, to an overall broadening of the Hubbard bands. In 
particular, this generates additional exponential tails in 
the vicinity of the band edges, whose extension is gov- 
erned by the variance a ~ \ZEpluq. In the adiabatic 
regime, where the boson frequency is small compared to 
the bandwidth (uio <§; D), and for large U, such tails are 
typically smaller than both the width of the Hubbard 
bands and the size of the gap (see figure l.a) 

The spectral function A(ek,oj), for U/D = 3 and 
a = 0.3Z?, is illustrated in figure 2. a at different values 
of £fc. This value of the variance corresponds, for exam- 
ple, to a moderate electron-boson coupling E P /D = 0.9, 
in the adiabatic regime luq = 0.1D and at zero temper- 
ature. The purely electronic case (a — 0) is also shown 
for comparison (dashed green lines). The electron-boson 
coupling strongly alters the lineshapes: well inside the 
Hubbard bands, there is a huge broadening of the peaks, 
which can be estimated in the large U limit (see ap- 
pendix) to Ar/r - 12a 2 /D 2 , and is of the order of 100% 
in the present example. However, the effect is even more 
dramatic in in the vicinity of the band edges. There, 
the scattering rate, which is roughly proportional to the 
spectral density, has a sharp (square root) dependence on 
the energy in the pure electronic case, causing a marked 
asymmetry of the peaks. The boson fluctuations con- 
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vert the sharp edge in a much smoother exponential tail, 
restoring a more symmetric lineshape for A(ef.,u>). Let 
us mention that the above-mentioned power law behavior 
at the edges of the Hubbard bands is not peculiar to the 
semi-circular DOS used in this example (see ref. [3ljp. 
Indeed, analogous results are obtained starting with a 
flat noninteracting DOS with a step-like edge, appropri- 
ate for two-dimensional systems. 

The position of the peaks is also affected by the 
electron-boson coupling, although to a weaker degree 
(this is partly in contrast with the results of refs. 
which only predict a broadening of the peaks, but no 
renormalization of the dispersion). In figure 2.b, we 
have reported the renormalized MDC dispersion Ef. . The 
main effect of the electron-boson coupling is an increase 
of the slope Vhe = dE^/deu; its relative variation is of 
the order a 2 /D 2 , as can be estimated from the large U 
expansion presented in the appendix, and is one order of 
magnitude smaller than the corresponding effect on the 
lincwidths. 



D. Isotope effect on the spectral properties 

We shall now analyze the consequences of a change 
in the boson frequency on the spectral properties. Un- 
less otherwise specified, we shall consider a relative shift 
Awo/wo = —6%, that can be achieved through the sub- 
stitution 16 — > ls O in the case of a lattice mode with 
predominantly oxygen character |32| . The value inferred 
from the photoemission experiments in ref. |17| on the 
mode at luq = 70meV is comparable with this value, 
|Acjo|/wo 7% , pointing to a strongly phononic char- 
acter of the bosonic excitation. 

By direct inspection of equation we immediately 
see that any isotopic effect will be strongly temperature 
dependent, on the scale of the boson frequency itself: 
at low temperatures, where a 2 = Epluq, a modification 
of u>o directly affects the distribution Q of the lattice 
displacements, and therefore modifies the spectral prop- 
erties described above. However, this effect is rapidly 
suppressed when the thermal fluctuations become dom- 
inant, in which case a 2 — 2EpT is independent on the 
boson frequency. In the following, we shall present the 
results at T = 0, where the IE is maximum. The IE at 
any temperature can be obtained straightforwardly by 
multiplying the results by an appropriate coefficient 



a(T) 



da 2 
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which is maximum at T — (a — 1), and rapidly drops 
at temperatures T > 0.2wo, where a ~ lu /3T. 

The inset of figure 2. a shows the IE on the peaks in 
A(ek,uj) at various close to the edge of the LHB. The 
reduction of the boson frequency leads to a weak shift 
of the peak position, corresponding to a slight reduction 
of the slope u^e, which is too weak to be observed at 
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FIG. 2: (color online) From top to bottom: a) energy scans 
of the spectral function A{ek,uj) in a Mott-Hubbard insu- 
lator with additional electron-boson interactions, at differ- 
ent tfe (equally spaced between the two edges = ±-D) 
(solid blue lines). The parameters are U/D — 3, Ep/D = 
0.9,uo/-D = 0.1, and energies are in units of D. The curves 
for ek/D = 1, 0.8, 0.6, 0.0 in the absence of electron-boson cou- 
pling are shown for comparison (dashed green lines). The in- 
set shows the effect of a shift Aloq/loq — —6% of the boson fre- 
quency on the spectral function at the same values of tk (the 
dotted red line is with the modified frequency), b) the dis- 
persion of the broad peaks deduced from momentum scans at 
constant energy, for U /D = 3 with (solid blue line) and with- 
out (dashed green line) electron-boson coupling. The slope 
of the noninteracting band is indicated. The inset shows the 
IE on the dispersion (the dotted red line is with the modified 
frequency) . c) Absolute value of the isotope effect on the scat- 
tering rate T — —ImT,, defined as AF = F[ojo + Ac^o] — r[wo], 
for different values of a 2 /D 2 = 0.06, 0.09, 0.12. 
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such moderate values of Ep (see also the inset of fig- 
ure 2.b, where we have reported the renormalized disper- 
sion close to the upper edge of the band). On the other 
hand, we know from the previous section that the effect 
of the electron-boson coupling is much stronger on the 
linewidths than on the dispersion, and indeed some re- 
duction of the peak widths is already visible in the plots 
of A(ek, oj). For a more quantitative analysis, we have re- 
ported in figure 2.c the IE on the scattering rate, defined 
as Ar = T[uj + AcJq] ^ r[fjj ]- It shows an interesting 
dependence on the energy: it is rather flat inside the Hub- 
bard bands, but exhibits pronounced peaks of width ~ a 
around the band edges, because this is the region where 
the spectra are mostly affected by the bosonic fluctua- 
tions. The variation of the linewidth in the flat region 
well inside the Hubbard bands is approximately given by 
(see appendix) 



Au D 



(11) 




and gives a direct measure of the strength of the electron- 
boson coupling. Note that although the precise numerical 
coefficient is specific to the semi-circular DOS considered 
here, analogous formulas can in principle be derived for 
any choice of the noninteracting band (in fact, a similar 
result also holds in the weak coupling regime 33J ) . 

Above the band edges, the IE decays exponentially, fol- 
lowing the fluctuation induced tails in the spectral den- 
sity. This gives rise to an extremely asymmetric peak 
of Ar at the band edge, which becomes more and more 
symmetric as a is increased (cf. the discussion on the 
spectral function in the previous section). For sufficiently 
large a > D, AT (a;) tends to a skewed gaussian, whose 
maximum is located in uj ~ —U/2 + 2<r, whose width is 
proportional to a and height scales as 



(12) 



III. POLARONIC SEMICONDUCTOR 

In this section, we use the approximate theory pre- 
sented above to address the spectral properties in a sys- 
tem with electron-boson interactions, but in the absence 
of electron-electron correlations. For this problem, exten- 
sive results are available in the litterature for the weak 
coupling regime (13 , and for the polaronic anti-adiabatic 
regime, where the boson frequency is assumed to be much 
larger than the noninteracting bandwidth 35]. We shall 
focus instead on the polaronic adiabatic regime (i.e. mod- 
erate to strong electron-boson couplings Ep > D and adi- 
abatic bosons ujq <C D), which is more often encountered 
in solids [3(| 153, HU , and for which a simple formulation 
of the spectral properties is not clearly established. 

For the problem of a single electron coupled to a disper- 
sionless boson, a complete characterization of the excita- 
tion spectrum has been given in references |l6l l34j , based 



on the Dynamical Mean Field Theory. In the adiabatic 
regime, the spectra are composed of one (or several, de- 
pending on the coupling strength) narrow features at low 
energy, equally spaced by ujq, coexisting with a continu- 
ous high energy background centered around the polaron 
binding energy Ep. As the coupling strength increases, 
the low-energy features are rapidly suppressed and the 
spectral weight becomes dominated by the high-energy 
incoherent background. It should be stressed that, con- 
trary to what happens in the anti-adiabatic limit, the 
high-energy features here are dispersive due to the strong 
hybridization with the free-electron states (see e.g. fig. 
14 in ref. [3). 

The high-energy incoherent excitations are well de- 
scribed by an equation analogous to , 



G(Q) = / dyP(y) 



1 



Go 1 (CI) 



(13) 



where the boson field obeys the same gaussian distribu- 
tion of eq. In fact, the above equation can be shown 
to be rigorously valid in the framework of the DMFT in 
the adiabatic limit ujo — [see ref. 0], eq. (46)]. A sim- 
ilar relation also holds for a system of spinless polarons 
at finite density leading to very similar results 

as in the single particle case presented here [ijj . 

Note that the chemical potential is undefined in the 
present single particle problem. To make direct contact 
with the results of the previous section, we shall interpret 
the results for the spectral function as corresponding to 
the lower "Holstein band" (i.e. to excitations at nega- 
tive binding energies as in a direct photoemission experi- 
ment), setting the chemical potential at the extremum of 
the polaron band, with the replacements u> = —CI + Eq 
and Cfc — > — efc in equation l|13|) . Here E < —D is 
the polaron binding energy taken from the DMFT so- 
lution, which tends to —Ep in the strong coupling limit 
E P /D-»1 (see fig. 4 in ref. [3). 

We shall first solve the coupled equations Ijl3(l and (J2J 
in a regime where the variance a of the boson field is 
smaller than the noninteracting bandwidth, which is typ- 
ically the case for moderate values of the electron-boson 
coupling. The opposite limit er > D, where the spectral 
density takes the form of a gaussian, multi-boson shakeoff 
peak, will be treated at the end of this section. 



A. Intermediate coupling regime 

The spectral density for Ep/D = 0.9 and wq/D = 0.1 
(a = 0.3D) is illustrated in figure l.b (blue dotted line). 
As in the correlated case, the boson fluctuations result in 
an overall broadening of the original band (green dashed 
line). The agreement with the DMFT result (black solid 
line) is excellent in the smooth region at high binding 
energies, while in the low energy region, the CPA clearly 
misses the detailed structure of the narrow peaks. How- 
ever, even there it gives a fair description of the inte- 
grated spectral weight, which is what one would mea- 
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FIG. 3: (color online) From top to bottom: a) Energy scans 
of the spectral function A{ej,,uj) in a polaronic semiconduc- 
tor, at different tu (equally spaced between the two edges 
€k = iD) within approximation eq. 1131 (solid blue lines). 
The parameters are Ep/D — 0.9,ujo/D — 0.1. The curves for 
et/D = 1,0.8,0.6, obtained using the exact DMFT solution, 
are shown for comparison (thin black lines). The inset shows 
the effect of a shift Aluq/ujo = —6% of the boson frequency 
on the spectral function at the same values of (dotted red 
line is with the modified frequency), b) the dispersion of the 
broad peaks deduced from momentum scans at constant en- 
ergy using DMFT (thin black line) and the approximate the- 
ory (solid blue line). The slope of the noninteracting band is 
indicated. The inset shows the IE on the dispersion obtained 
in DMFT (the red dotted line is with the modified frequency), 
c) Absolute value of the isotope effect on the scattering rate 
Ar using DMFT (thin black line) and the approximate the- 
ory (solid blue line) . Data from exact DMFT are shown after 
convolution with a gaussian filter of width 0.05D. 



sure experimentally in the presence of a sufficient energy 
broadening. 



1. High-energy features 

The results for the spectral function are illustrated in 
figure 3. a. As in the correlated case treated in the previ- 
ous section, the high energy part of the spectra is char- 
acterized by broad features, whose dispersion (fig. 3.b) 
sensibly deviates from the noninteracting case. The vari- 
ation of the scattering rate under a shift of the boson 
frequency (fig. 3.c) is also very similar to the correlated 
case, corroborating the fact that the magnetic fluctu- 
ations only play an indirect role in the isotope effect: 
Ar(w) is rather flat at the center of the band, and at- 
tains its maximum in a narrow region of width a around 
the band edges. 

For a more quantitative understanding, analytical ex- 
pressions can be obtained by an expansion to lowest or- 
der in the variance, as was done previously for the Mott- 
Hubbard insulator. In the present case the coupled equa- 
tions l|13|) and @ become 



S 
G 



<j G 



1 



n-{D 2 /A + a 2 )G 



(14) 
(15) 



which is valid inside the band, far from the edges. The 
Green's function has the noninteracting form, but with 
a renormalized bandwidth given by D -> VD 2 +4a 2 . 
The slope of the dispersion of the high-energy features 
is renormalized accordingly: v/ le = dE^/de^ = 1/(1 — 
2<r 2 /D 2 ), i.e. it does not coincide with the value in the 
absence of interactions, as can be seen in figure 3.b. The 
boson induced variation is of the same order as what was 
calculated in section 11. C in the presence of electronic 
correlations. 

From the first equation we see that the scattering rate 
within the band is proportional to the spectral density, 
T(f2) ~ Tia 2 N*(Q). Note that it is weaker than in the 
correlated case, due to the absence of magnetic disorder. 
Nevertheless, the modification of the scattering rate un- 
der a shift of the boson frequency is of the same order as 
in the previous case [although with a smaller prefactor, 
cf. eq. ijnj], namely 



Ar 



2E F 



(16) 



being directly proportional to the strength of the 
electron-boson coupling. 

The existence of a maximum of Ar near the band edge 
(figure 3.c) also compares well with the DMFT result, 
although its actual position is slightly shifted to higher 
binding energies. This can be understood by observing 
that the true edge of the incoherent dispersion is at lu = 
—loq (not at u> = 0), which marks the boundary between 
the high-energy and low-energy regions in the excitation 
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spectra (see below). Note that Ar rapidly drops in the 
region — lo q < u) < 0, where Y itself is extremely small 
due to our assumption of dispersionless (gapped) bosons. 



2. Low-energy features 

The low-energy part of the spectral function is shown 
in the inset of figure 3. a. The DMFT result shows that 
the broad incoherent peak progressively disappears when 
the band edge is approached (i.e. at low momentum 
transfers), while a narrow "quasi-particle" peak arises at 
binding energies \u)\ < luq. The evolution of such peak- 
hump structure, which is characteristic of the intermedi- 
ate coupling regime (in the strong coupling regime, the 
narrow features are too weak to be observed) causes a dis- 
continuous jump, or kink, in the dispersion, which clearly 
separates the high and low-energy regions with different 
slopes, as illustrated in figure 3.b (see also figure 3 in 
reference The isotope effect on the kink region is 

shown in the inset. Note that in this plot, the isotope 
shift vanishes at to — by definition, since the origin of 
energies has been shifted to coincide with the band edge 
(see the discussion at the beginning of this section). 



B. Strong coupling limit 



At the center of the polaron peak, the self-energy 
E(w) = -u)-E P - 1/G(u) tends to 

Z(u) = -(uj + E P )(l-~\-i^- (19) 

As a consequence, the dispersion of the broad peaks in 
A(ek, uS) tends to = —Ep + (tt/2) e^, which defines an 
apparent bandwidth = -kD for the incoherent features, 
sensibly larger than the noninteracting value 2D [this 
should not be confused with the width of the polaron 
peak in the momentum-integrated spectral density of eq. 
(|17f) . which is governed by the variance a]. Correspond- 
ingly, the slope of the high-energy dispersion saturates 
to a finite value in the strong electron-boson coupling 
limit, which is independent of the coupling strength, and 
is larger than the noninteracting value (cf. the situation 
in the Mott- Hubbard insulator at strong U, section II. B). 

In the energy interval spanned by the dispersion Ek, 
the scattering rate is roughly constant and directly pro- 
portional to a. Its variation under a change of the boson 
frequency is given by 

Ar n^r 

which is similar to the behavior encountered near the 
edges of the Hubbard bands [cf. eq. (12)]. 



At extremely large values of the coupling strength 
Ep/D > D/ujq ^> 1 (or at sufficiently high tempera- 
tures T I D > D/Ep), the variance of the boson field can 
become comparable with the noninteracting bandwidth 
{cr~>D). The latter can therefore be neglected in eq. 
replacing Gq 1 = fl. Note that this does not correspond 
to the usual atomic, anti-adiabatic limit [351 ] . where it is 
assumed from the beginning that D — > is the smallest 
energy scale in the problem, resulting in dispersionless 
high energy features. The present theory is valid in the 
opposite limit, D > uq, which is more often realized in 
solids (3(J 153, HH]. Due to the large transfer integrals 
between molecular units, the discrete shakeoff spectrum 
characteristic of isolated molecules is converted here into 
a continuous gaussian spectral density 0] , and a sizeable 
high-energy dispersion is recovered. 

We recognize from eq. (113(1 that the spectral density in 
this case coincides with the gaussian distribution itself, 



N*(lo) =P{uj + E p ) 



(17) 



whose width is governed by the variance a (we have 
performed the shift VI = —to — Ep using the fact that 
Eq — > — Ep in the strong coupling limit). The full 
Green's function G can be read directly from eq. 10: 



Ep + iS 



(18) 



IV. CONCLUDING REMARKS 

In this work, we have presented an approximate ana- 
lytical theory which addresses the high energy spectral 
properties in systems characterized by a strong electron- 
boson coupling, both in the presence and in the absence 
of electronic correlations. Concerning the electron-boson 
interaction alone, the present approach gives accurate re- 
sults in the adiabatic regime (i.e. opposite to the stan- 
dard strong-coupling polaron theories), where its validity 
can be controlled by direct comparison with the results 
of the Dynamical Mean Field Theory. For the correlation 
part, on the other hand, it reduces to the CPA treatment 
of reference |19| , which qualitatively accounts for the in- 
coherent high-energy excitations located in the upper and 
lower Hubbard bands. 

Although the microscopic mechanisms in the cuprates 
certainly go beyond the simple model and approxima- 
tions presented here, our results reproduce at least quali- 
tatively several characteristics of the observed photoemis- 
sion spectra, such as the existence of broad peaks with 
a sizeable momentum dispersion, which sharpen and be- 
come more asymmetric (in the EDC scans) as the band 
edge is approached. 

More specifically, we have calculated the isotope ef- 
fect on the high-energy spectral features. In the present 
framework, where the bosonic and magnetic fluctuations 
are effectively decoupled, the overall behavior in the Mott 
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insulator and in the polaronic semiconductor is qualita- 
tively similar. The most notable result is the existence 
of a strongly energy dependent IE on the linewidths, 
which is maximum in a narrow energy interval of width 
a ~ \/Epu;q near the band edges, because this is the 
region where the excitation spectrum is mostly affected 
by the electron-boson interaction (new spectral weight is 
created there, due to the presence of boson fluctuations) . 

The present approach also predicts a strong temper- 
ature dependence of the IE, which should be strongly 
suppressed when the temperature reaches some fraction 
of the boson frequency (typically T > 0.2wo). This oc- 
curs because at high temperatures, the fluctuations of 
the boson field are dominated by thermal effects, which 
are independent on the boson frequency, and is in no 
way related to the existence of a temperature depen- 
dent electron-boson coupling. Note that, according to 
the above general arguments, an analogous energy and 
temperature dependence of the IE can also be expected 
in more accurate treatments of the electron correlations. 

The results of section II can be tentatively compared 
with the experimental results of references [IT], 13 ■ If we 
associate the variance a with the observed width ~ 0.2eV 
of the active IE region, and take the value loo = 0.07eV 
for the boson energy, a "polaron" binding energy Ep ~ 
O.beV is obtained. Assuming a noninteracting bandwidth 
of the order of leV places this value in the intermediate 
electron-boson coupling regime. The same value of a is 
also compatible with the magnitude of the observed IE on 
the scattering rate. From the frequency softening Aloq = 
5—10meV deduced from the shift of the kink energy, the 
theory predicts a decrease of the scattering rate around 
the band edges of the order Ar ~ aAujo/uj ~ 10 — 
30meV (see the end of section II. D), in agreement with 
the experimental observations of ref. [l8| . 

It should be stressed that the present theory, based on 
a momentum independent electronic self-energy, clearly 
fails in addressing the strongly anisotropic dispersion ob- 
served in the cuprate superconductors, which demands to 
go beyond a local (and classical) treatment of the bosonic 
and magnetic fluctuations. In particular, the edges of 
the Hubbard bands at low temperature are strongly af- 
fected by the "spin density wave" dispersion, especially 
at intermediate values of U. Another point which is be- 
yond the range of validity of the present approach is the 
enhancement of the low energy IE due to proximity to 
a Mott metal-insulator transition. In this case, even a 
weak electron-boson interaction may produce a huge ef- 
fect when the boson frequency is varied |15(. 
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APPENDIX: IE IN THE LARGE U LIMIT. 



the UHB in equation J2J for sufficiently large U. In the 
case of a semi-circular DOS, we obtain 



G 



1/2 



lo - D 2 G/A + U/2 

1/2 

lo- {D 2 /A + 2a 2 )G + U/2 



[lo - D 2 G/A + U/2] 



(21) 



where we have used the fact that lo — D 2 G/4 + U/2 = 
1/(2G) to lowest order in a 2 . The solution for G reads 



G(co) = -L 



U/2 -yf{w + U/2) 



2 £)* 2 



(22) 



with D* 2 = D 2 /2 + 4cr 2 . From the self-consistency rela- 
tion G _1 =lo- D 2 G/4 - Ewe obtain 



E(w) ~ -lo - U + (D 2 /A + 4cr 2 ) G{lo) 



(23) 



The scattering rate T(co) = — /mE(u) is therefore di- 
rectly proportional to the spectral density at this energy. 
At the center of the band, it is given by 



Dy/2 



Its isotope effect at T = 



12a 2 



D 2 



Alo q = '^-D 



(24) 



(25) 



gives a direct measure of the strength of the electron- 
boson interaction. 

b. Fluctuation induced tails. In the gap region be- 
tween the Hubbard bands, additional spectral weight is 
created by the boson fluctuations. For small er and suf- 
ficiently far from the band edges, both terms in can 
be replaced by their atomic counterparts leading to the 
following exponential decay (for the LHB) 



7T \LO 



U 2 /4} 2 



8a 2 



(26) 



(note that this formula is not valid close to the band edge, 
i.e. where the IE is maximum). 

A simpler result is obtained for sufficiently large a, i.e. 
in the strong electron-boson coupling regime U 3> cr 2 ^ 
D. In this case, the variation AT(lo) of the scattering rate 
under a shift of the boson frequency takes the form of a 
skewed gaussian, which is maximum at lo ~ — U /2 + 2a, 
where 



(27) 



a. Inside the Hubbard bands. Well inside the LHB, 
we can expand © for small a 2 , and neglect the effect of 
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